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I. INTRODUCTION 


Due to the aging of buildings or even the change in use 
of the structure, reinforcing RC structures can become the 
most economical and viable option to be carried out. 
According to Souza and Ripper [1], the reinforcement of 
columns is always more problematic than beams and slabs. 
The columns are the last elements of support of the 
structure before the loads reach the foundations and have 
to obey safety and economic criteria in any situation. 


Among the existing methods, the reinforcement of 
columns with CFRP (Polymer Reinforced with Carbon 


Fiber) has become increasingly common. Saadatmanesh 
and Ehsani [2] state that the growth in the use of 
composites is mainly due to their high tensile strength, 
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Abstract—Carbon Fiber Reinforced Polymer (CFRP) reinforcement of 
reinforced concrete (RC) columns has become increasingly common. This 
work presents a finite element computational model capable of 
representing the structural behavior of CFRP-confined square RC 
columns. Data and results from experimental research available in the 
literature are used to assess the ability of the proposed model to represent 
the test results. It is possible to observe that the models were able to 
represent, with adequate precision, the structural behavior of columns 
with different configurations of longitudinal reinforcement, stirrups 
spacing, number of CFRP layers, and type of load. 


high modulus of elasticity, and excellent corrosion 
resistance. 


This reinforcement technique in columns works due to 
the lateral confinement promoted by the CFRP. However, 
according to Carrazedo [3], confinement is only 
appropriately enabled with the development of 
microcracks at a high-stress level since the lateral 
deformation of concrete is very small before the beginning 
of microcracking. Therefore, this reinforcement technique 
is classified as a passive form of confinement. 


Several experimental works in the literature address the 
reinforcement of RC columns with CFRP wrap. However, 
elaborating all the procedures of experimental research can 
become onerous. Therefore, using computational methods 
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capable of representing the structural behavior of columns 
can be considered an interesting alternative. 


Computer simulations of concrete structures, for the 
most part, are with the finite element method (FEM). Real 
et al. [4] explain that in the FEM, the structure is divided 
into small elements that are linked together through nodal 
points. This method solution is determined through the 
displacements, which are the unknowns of the problem. 
The internal stresses of the element and the reactions at the 
supports are calculated from these displacements. 


In Brazil, the Brazilian Standard (NBR) that deals with 
reinforced concrete elements’ design and structural 
behavior is the NBR 6118 [5]. However, this standard still 
does not provide prescriptions regarding the design and 
construction procedures of CFRP _ reinforcements. 
Therefore, conducting research in this area is essential in 
the national context. 


Most Brazilian research on carbon fiber reinforced 
structures is based on international normative 
recommendations. Therefore, this research aims to 
demonstrate the computational simulation of RC columns 
reinforced with CFRP to contribute to the scientific 
community as a technical asset for future Brazilian 
standardization of this reinforcement technique. 


Il. METHOD 


This work was carried out using the experimental 
results of research by Wang et al. [6], Atheer et al. [7], and 
Souza [8] as a reference in terms of comparison to validate 
the computational model used. The computational 
modeling was performed using the Ansys Workbench 17.1 
software, which performs simulations using the finite 
element method. 


2.1 — Finite element analysis (FEA) 


The SOLID65 element was used to represent the 
concrete. According to ANSYS [9], SOLID65 is a 3D 
element, with eight nodes, with a material model capable 
of cracking under tension and failure by crushing in 
compression and having the ability to suffer creeping. This 
element supports up to four different materials for each 
finite element, with one being the principal material and 
the reinforcement bars distributed inside the element 
volume, being ideal for representing the reinforced 
concrete 


The linear element LINK180 represented the 
longitudinal reinforcement and stirrups of the columns. 
According to ANSYS [9], the LINK180 element can be 
used to model cables and bars and has two nodes with 
three degrees of freedom at each node. This element resists 
only axial forces and considers the steel reinforcement's 
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elasticity, plasticity, and creep characteristics. In addition, 
its cross-section can increase or decrease when the element 
is subjected to compression or tension. 


The SHELL 181 element was used for the CFRP, 
which is capable of considering the change in the thickness 
of the CFRP in nonlinear analysis, a characteristic that 
makes it able to evaluate the internal stress distribution 
during the study of the structure where it is associated. 
According to ANSYS [9], it is a four-node element with 
six degrees of freedom at each node (translations in the x, 
y, and z directions and rotations about the x, y, and z axes). 
Fig. 1 represents the elements used. 


‘STEEL; LINKISO 


CONCRETE: SOLIDE6S 


$ 
CFRP: SHELL181 


Fig. 1- Types of the elements used in the FEA 


The failure mode of the SOLID65 element used is 
based on the criterion of William and Warnke [10]. These 
equations consider both the crushing and cracking failure 
of the concrete. Poisson's coefficient was 0.2 in all 
columns, while the hydrostatic, shear transfer coefficients 
for closed cracks and shear transfer coefficients for open 
cracks were defined as 1.0, 0.3, and 1.0, respectively. 


A discrete model was used to idealize the stiffness of 
the steel. According to Ayala [11], it is a model that 
requires that the nodes of the SOLID65 and LINK180 
elements coincide. This element combination is the ideal 
model for the constitutive equations of William and 
Warnke [10]. 


A mesh convergence test was performed to determine 
the best mesh configuration. For the columns of Wang et 
al. [6], meshes of 7,625 mm were used in x and y 
directions for concrete, while for CFRP, a mesh of 20 mm 
was employed. In the columns of Atheer et al. [7], meshes 
of 7.5 mm were adopted in x and y directions for concrete, 
and a mesh of 12.5 mm was defined for CFRP. Finally, in 
the columns of Souza [8], were used meshes of 6.25 mm in 
the x-direction and 28.52 mm in the y-direction for 
concrete and a mesh of 12.5 mm was defined for CFRP. 
The meshes used are shown in Fig. 2 - 4. 
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CONCRETE CRP 


Fig. 2 — Mesh of the columns by Wang et al. [6] 


CONCRETE CFRP 


Fig. 3 — Mesh of the columns by Atheer et al. [7] 


CONCRETE ERP 


Fig. 4 — Mesh of the columns by Souza [8] 
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The analysis was performed considering the lower face 
of the columns fixed in the Y direction to represent the 
models’ reaction. In contrast, the upper face was 
considered free to move. A displacement in the Y direction 
was applied in this region to represent the load action. This 
action was subdivided into "substeps" adjusted for each 
model. Like the mesh convergence test, the appropriate 
number of "substeps" depends on tests and the geometry of 
each model. 


The Newton-Raphson method, which was used in this 
analysis, is one of the ways used by ANSYS [9] to solve 
nonlinear analyses. According to Pivatto [12], this method 
approaches the equilibrium trajectory of the structure by its 
tangents until reaching convergence. The displacement 
increments were applied until the column load capacity 
had been exceeded. 


IW. RESULTS AND DISCUSSION 
3.1 — Experiments by Wang et al. [6] 


The experimental work by Wang et al. was carried out 
to study the influence of the variation of the stirrups 
spacing and the number of layers CFRP in columns of 
larger dimensions than are generally found in the 
literature. The authors tested 34 columns with two 
identical specimens; there were 17 column configurations. 
One column of each format was tested for a monotonic 
load, and the other was tested for a cyclic load. Fig. 5 
presents the dimensions of the columns and the details of 
the longitudinal reinforcement and stirrups. 


A i 
$6@120 
26 @80 | ($6@60) 
($6@40) a 
od 
B B 
66@30 + 6@30 ! 
305 204 


Fig. 5 — Specimen dimensions and reinforcement 
details of the columns by Wang et al. [6] 
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Table 1 — Characteristics of the materials used by Wang et al. [6] 
Material Diameter/thickness (mm) Yield/peak strength Ultimate strength Elastic Modulus 
Longitudinal steel 10 312 484 200,000 
12 340 520 200,000 
Stirrupssteel 6 397 623 200,000 
CFRP wrap 0.167 - 4,340 240,000 


The models are characterized as follows: The letters S1 
and S2 represent the dimension of the analyzed cross- 
section, and the letter H represents the stirrup rate used, 
which can be 0 for simple concrete, 1 or 2 for a stirrup 
steel rate of 0.5% or 1% respectively. The letter L 
represents the number of CFRP layers, which can vary 
from 0 to 3 layers, and finally, the letter M represents the 
type of load used, which would be monotonic. Fig. 5 
illustrates the dimensions of the tested columns and the 
reinforcement details. For example, the model S1H1L3M 
has a square cross-section of 305 mm, a stirrup steel rate of 
0.5%, 3 layers of CFRP, and has been subjected to a 
monotonic load. 


The columns received a rounding of the corners to 
mitigate the effects of stress concentration at the corners. It 
was used a rounding with a radius of 30 mm for the S1 
specimens and 20 mm for the S2 specimens, keeping a 
ratio of twice the rounding radius per cross-section side at 
0.2 so that the two models receive the same influence. The 
compressive strength of the concrete was 25.5 MPa on the 
day of the tests. 


The stress was obtained by dividing the maximum load 
by the cross-sectional area as in the experimental results. 
The deformation was calculated by the ratio between the 
axial displacement and the height of the analyzed column. 
For the simulation of the columns by Wang et al., axial 
displacements on the Y-axis were applied to the upper face 
of the columns. This displacement was subdivided into 
200 substeps. 


In terms of comparison, only the results of simulations 
performed with monotonic loads, i.e., continuous loading, 
were used. On the other hand, cyclic loads are variable 
loads, i.e., successive loading and unloading until the 
element under study fails. Furthermore, the authors did not 
find significant differences between the results obtained 
between the different types of loading. Table 2 
demonstrates the comparison between experimental and 
computational results. 
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The authors reported that the confinement distribution 
of the CFRP was not uniform around the perimeter of the 
columns. The CFRP was not efficient in confining the 
concrete in the region in the middle of the lateral face of 
the columns. On the contrary, the corners were responsible 
for promoting the observed confinement. Therefore, the 
mean strain value obtained overestimated the failure stress 
and the effectiveness of the confinement of the CFRP. 
Wang et al. [6] concluded that the corners were the most 
critical points and that the CFRP strain efficiency factor 
was about 0.4 for the larger columns (S1) and 0.6 for the 
smaller columns (S2). 


As in the experiments by Wang et al. [6], the computer 
models showed a more significant deformation in the 
middle of the lateral faces of the CFRP. Furthermore, the 
models with a smaller cross-section presented a greater 
lateral deformation when compared to the models with a 
larger cross-section, which is in line with the results given 
by Wang et al. [6]. Figs. 6 — 8 represent the strain on 
CFRP of the larger (S1) and smaller (S2) section models. 


The other columns showed the same strain distribution 
pattern in Figures 6 - 8. Wang et al. [6] reported that the 
models suddenly failed in the mid-height region regarding 
the failure mode. Furthermore, the RC columns also 
showed buckling in the longitudinal reinforcement bars 
between the stirrups, and the concrete had been severely 
crushed. Fig. 9 represents the stress in the longitudinal 
steels and stirrups of the computational models S1H1L3M 
and S2H2L2M. 


0,0029167 Min 


(a) 
Fig. 6 — Lateral strain SIHOL2M (a) 1° layer (b) 2° layer 
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Table 2 — Comparison between the results of Wang et al. [6] and the proposed model 


Specimen 


S1H1LOM 
S1H2LO0M 
S1HOLIM 
S1HOL2M 
S1H1ILIM 
S1H1L2M 
S1H1L3M 
S1H2L2M 
S1H2L3M 
S2H1LOM 
S2H2L0M 
S2HOLIM 
S2HOL2M 
S2H1LIM 
S2H1L2M 
S2H2L1M 
S2H2L2M 


00038201 Min 


(a) 


Fig. 8 — Lateral strain SIHIL3M (a) 1°layer (b) 2° layer 
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Experimental 


Axial strain 


(mm/mm) 
0.00373 
0.00412 
0.00313 
0.00393 
0.00428 
0.00434 
0.00428 
0.00416 
0.00529 
0.00364 
0.00406 
0.00391 
0.02287 
0.00704 
0.03588 
0.00465 
0.04230 


QONA 
000099871 Min 


(a) 


Axial stress 


(MPa) 
32.10 
34.70 
29.40 
32.30 
35.10 
34.90 
36.90 
35.50 
37.20 
29.90 
32.50 
28.70 
31.40 
35.50 
40.00 
33.00 
40.80 


pp 0088207 Mex 
FB scone 
gone 
0005055 
D00325 
noo 2935 
00011574 
Damin 
¢.00097536 Min 


(b) 


Fig. 7 — Lateral strain S1HIL3M (a) 1° layer (b) 2° layer (c) 3° layer 


Computational 


Axial strain 


(mm/mm) 
0.00241 
0.00249 
0.01090 
0.01311 
0.00459 
0.00361 
0.00448 
0.00426 
0.00535 
0.00359 
0.00327 
0.00817 
0.01307 
0.00653 
0.01471 
0.00474 
0.01634 


000093118 Min 


(c) 


Axial stress 


(MPa) 
32.79 
34.77 
28.60 
29.73 
36.43 
36.70 
37.60 
37.50 
38.50 
30.24 
32.89 
28.69 
31.98 
35.00 
39.996 
34.77 
40.76 


Variation 
between 


stress 


+2.15% 
+0.20% 
-2.712% 
-7.96% 
+3.79% 
+5.16% 
+1.90% 
+5.60% 
+3.50% 
+1.14% 
+1.20% 
-0.25% 
+1.85% 
-1.41% 
-0.01% 
+5.36% 
-0.10% 


Fig. 9- Stress in the steels (a) SIHIL3M(b) S2H2L2M 


Page | 75 


Silva et al. 


As seen in Fig. 9, the longitudinal steels reached the 
rupture stress specified in Table 1, presenting the same 
failure characteristics as the RC columns tested by Wang 
et al. [6]. In addition, the authors reported that some of the 
columns with the largest section (S1) showed post-peak 
softening, which can also be observed in models 


S1HOL1M 

50 
© 40 
= 
= 30 
wn 
g 
*% 20 
& 
ž 10 

(0) 

0.00 0.01 0.02 0.03 0.04 0.05 
Axial strain (mm/mm) 
Experimental Computacional 
S1H1L1M 

50 
© 40 
= 
= 30 
wn 
g 
tn 20 
Xs 
3 10 

(0) 

0.00 0.01 0.02 0.03 0.04 0.05 


Axial strain (mm/mm) 


Experimental Computational 


S1H1L3M 


fo 


0.00 0.01 0.02 0.03 


N w p 
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S1H1L1M and S1H1L2M in the stress-strain curves of 
Fig. 10. In addition, the smaller section models (S2) 
showed an upward behavior, as can be observed in the 
S2H1IL2M and S2H2L2M models, which was also 
reported by Wang et al. [6] in their experiment. 


S1HOL2M 
50 
= 40 
= 
J 30 
n 
g 
% 20 
S 
%ž 10 
0 
(0) 0.01 0.02 0.03 0.04 0.05 
Axial strain (mm/mm) 
Experimental Computational 
S1H1L2M 
50 
© 40 
= 
= 30 
nn 
2 
ï 20 
S 
3 10 
(0 
(0) 0.01 0.02 0.03 0.04 0.05 
Axial strain (mm/mm) 
Experimental Computational 
S1H2L2M 
50 
© 40 
2 
z 30 
n 
g 
t 20 
T 
ž 10 
0 
0 0.01 0.02 0.03 0.04 0.05 


Axial strain (mm/mm) 


Experimental 


Computational 
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S1H2L3M 
50 
© 40 
= 
J 30 
n 
2 
Į 20 
is) 
3 10 
(0) 
0.00 0.01 0.02 0.03 0.04 
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Experimental Computational 
S2HOL2M 
50 
© 40 
2 
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L 
g 20 
T 
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= 
z 30 
wn 
g 
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S 
ž 10 
0 
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S2HOL1M 
50 
40 
30 
20 
10 
(0) 
0.00 0.01 0.02 0.03 0.04 0.05 
Axial strain (mm/mm) 
Experimental Computational 
S2H1L2M 
50 
40 
30 
20 
10 
0 
(0) 0.01 0.02 0.03 0.04 0.05 
Axial strain (mm/mm) 
Experimental Computational 
S2H2L2M 
50 
40 
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20 
10 
o 
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Experimental Computational 


Fig. 10 — Stress-strain graphs of the computational model for the columns of Wang et al. [6] 
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3.2 — Experiments by Atheer et al. [7] 
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Table 3 — Characteristics of the columns by Atheer et al. [7] 


Specimen Section Height Longitudinal Stirrups(mm) N° CFRP Eccentricity 
(mm) (mm) steel (mm) wrap (mm) 
SFO 150 x 150 800 4@ 12 Ø6 c/50 2 0 
SF20 150 x 150 800 4@ 12 Ø6 c/50 2 20 
SF40 150 x 150 800 4@ 12 Ø6 c/50 2 40 
The authors carried out an experimental study (a) 
proposing a new form of reinforcement. For this, 


experimental simulations of CFRP-confined square RC 
columns were carried out to serve as a benchmark. 
Subsequently, 
reactive powder was performed in models identical to the 


circularization by concrete based on 


square ones, increasing the cross-section and changing the 
shape of the section of the column studied. 


For the simulations, only the reference specimens were 
selected, that is, the CFRP-confined square models without 
performing the circularization, which were subjected to 
concentric and eccentric loads. In addition, all the columns 
have the same configuration, differing only by the 
eccentricity used. The characteristics of these columns can 
be seen in Fig.11. 


4N12 mm 
R6@50 mm 


Fig. 11 — Cross-section of columns (a) S and (b) SF 
(Atheer et al. [7]) 


For the columns by Atheer et al. [7], the displacement 
applied to the upper face was subdivided into 75 substeps 
for model convergence. The corners received a rounding 
radius of 20 mm to alleviate the stress concentration in this 
region. Table 3 presents the details of the columns tested 
by Atheer et al. [7], while Table 4 presents the 
characteristics of the materials used. 


Table 4 — Characteristics of the materials used by Atheer et al. [7] 


Material Average Ultimate strength Elastic Average Ultimate 
compressive (MPa) Modulus tensile (MPa) strain (%) 
(GPa) 
(MPa) 

Concrete 25 - - - - 

Longitudinal steel - 634 - 2 z 

Stirrups steel - 613 - - - 

CFRP wrap - - 291 2,668 2.47 


It was reported in the research that the CFRP-confined 
model subjected to concentric load (SFO) failed due to 
fragmentation of the cover due to stress concentration near 
the corners in the mid-height region. Fig. 12(a) shows 
these stress concentrations. In addition, it was reported that 
the longitudinal reinforcement bars buckled, as shown in 
Fig. 12(b). 


In the specimen subjected to an eccentric load of 20 


mm (SF20), fragmentation of the cover on the 
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compression side was observed due to stress concentration. 
The buckling of the reinforcement bars on the compression 
side of the column was also observed. Figs. 13 (a) and (b) 
illustrate the stresses on concrete and steel reinforcements. 
In addition, there were cracks between the CFRP strips on 
the tensioned side with CFRP rupture on the compression 
side. 


The failure mode of the model subjected to 40 mm 
eccentricity (SF40) was similar to that of the SF20 model, 
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showing fragmentation of the concrete on the compression 
side and cracks between the CFRP strips on the tensioned 
side. Fig. 14 illustrates the stresses in concrete and 


reinforcement steel bars. 


The summary of the simulation results can be seen in 
Table 5, where the maximum load of the computational 
models is compared to the results obtained by Atheer et al. 


[7]. 


-75,922 Min 


(a) 


ONNRERBEL £2 


Fig. 12 — SFO (a) stress in the concrete (b) stress in the 


steel reinforcements 


-45,014 Min 


(a) 


Fig. 13 — SF20 (a) stress in the concrete (b) stress in 
the steel reinforcements 


1,2306 Max 
-2,4118 
-60541 
-9,6965 
-13,339 
-16,981 
-20,624 
-24,266 
-27,908 
-314,551 Min 


(a) 


= 
= 
= 
= 
==] 
SS 
= 
i 


Fig. 14 — SF40 (a) stress in the concrete (b) stress in 
the steel reinforcements 


Table 5 - Comparison between the results of Atheer et al. [7] and the proposed model 


Specimen Experimental 
Axial 
displacement 
(mm) 
SFO 30.50 
SF20 9.10 
SF40 8.9 


www.ijaers.com 


Axial 
force 
(KN) 
1,462 
650 
429 


Computational 
Axial Axial 
displacement force 
(mm) (KN) 
30.00 1,466.2 
3.03 670.11 
2.70 462.69 


Variation 
between axial 


forces 


+0.29% 
+3.09% 
+7.85 % 
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3.3 — Experiments by Souza [8] 


Table 6 - Characteristics of Souza's [8] columns 


Specimen Section Height Longitudinal steel Average N° 
(mm) (mm) (mm) compressive CFRP 
(MPa) wrap 
“GI = 150x150 7B 400o O a2 o 
G2 150 x 150 713 4Ø10 23.0 - 
G3 150 x 150 713 40 12.5+4@16 41.2 - 
GIR 150 x 150 713 4@ 10 41.2 1 
G2R 150 x 150 713 4@ 10 23.0 1 
G3R 150 x 150 713 40 12.5+4@16 23.0 1 
In his doctoral thesis, Souza [8] carried out designed one. This situation can be found in structures of 
experimental tests of 6 columns, being 3 CFRP-confined. very advanced ages. Table 7 shows the characteristics of 
The columns had the same geometric characteristics, the other materials used in Souza's experiment [8]. 


varying only the longitudinal reinforcement steel rate and For the columns by Souza [8], the displacement applied 


the compressive strength of the concrete. The to the upper face was subdivided into 50 substeps for 


characteristics of the columns can be seen in Table 6. model convergence. The corners received a rounding 
In the G3R model, the author used concrete with lower radius of 25 mm to alleviate the stress concentration in this 
strength than his reference model to simulate a situation in region. 


which the column had a bearing capacity lower than the 


Table 7 - Characteristics of the materials used by Souza [8] 


Material Ultimate strength Elastic Modulus Average tensile Ultimate 
(MPa) (GPa) (MPa) strain (%) 
Longitudinal steel 615 - - - 
Hoop steel 797 - - - 
CFRP wrap - 231 4,205 1.8 
It was reported by Souza [8] that the failure mode of stresses (Fig. 16) was also shown to be similar to the 
the unreinforced models occurred initially in the cover reports by Souza [8]. Therefore, the computational model 
region. Later there was the crushing of the concrete core, failure also occurred due to the crushing of the concrete 
with buckling of the longitudinal steel bars. The models core of the column. Table 8 combines the experimental 
CFRP-confined showed rupture of the concrete core but and computational simulation results, comparing the 
without showing fracture of the CFRP. results obtained. 


The computer models showed similar behavior. The 
CFRP did not reach its maximum tensile capacity in any of 
the models (Fig. 15). The distribution of compressive 
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67 
177,17 

| 147,64 
118,11 
88,586 
59,057 
29,529 
-1,0806e-11 Min 


-1,5409e-11 Min 


-1,8534e-11 Min 


(c) 


(a) 


Fig. 15 — Stress in the CFRP (a) GIR (b) G2R (c) G3R 


-34,411 Min 


Fig. 16 — Stress in the concrete (a) GIR (b) G2R (c) G3R 


Table § - Comparison between Souza's [8] results and the proposed model computational 


Specimen Experimental Computational Variation 
~~ Axial strain Axial Force Axial strain. Axial force between axial 
(mm/mm) (KN) (mm/mm) (KN) forces 
Gl 0.002423 1013 0.002104 1042.3 + 2.89% 
G2 0.001852 613 0.002524 600.4 - 2.05% 
G3 0.002061 1444 0.002384 1377.9 - 4.58% 
GIR 0.004600 1417 0.008415 1324.2 - 6.55% 
G2R 0.005586 927 0.007010 861.26 - 7.09% 
G3R 0.005653 1303 0.004909 1401.3 + 7.54% 
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IV. CONCLUSION 


The proposed computational model demonstrated 
satisfactory performance. It showed good behavior in the 
stress-strain graphs compared to the columns of Wang et 
al. [6]. Still, on the columns of Wang et al. [6], it 
presented good agreement in the ascending branch of the 
models of smaller cross-sections and the post-peak 
softening of some models of the larger section. 


In the Atheer et al. [7] columns, the simulation 
presented maximum loads close to the experimental 
results, mainly in the concentric load model. However, it 
showed a different axial displacement in the models with 
eccentric loads. 


Finally, in the columns tested by Souza [8], the model 
also had satisfactory results, demonstrating that the high 
rate of longitudinal steel of the G3R column reduced the 
lateral expansion of the column, which reduced the 
contribution of the CFRP in the final axial capacity of the 
column. 
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